Strong-weak coupling self-duality in the two-dimensional quantum phase transition of 

p + ip superconducting arrays 
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The 2D quantum phase transition that occurs in a square lattice of Josephson-coupled p ± ip 
superconductors is an example of how four-body interactions in d = 2 reproduce nonperturbative 
effects caused by two-body interactions in d = 1. The ordered phase has unconventional "bond 
order" of the local T-breaking variable. This problem can be analyzed using an exact self-duality; 
this duality in classical notation is the 3D generalization of the Kramers- Wannier duality of the 2D 
Ising model, and there are similar exact dualities in dimensions d > 3. We discuss the excitation 
spectrum and experimental signatures of the ordered and disordered phases, and the relationship 
between our model and previously studied behavior of 2D boson models with four-boson interactions. 

PACS numbers: 74.50-(-r 74.72-h 



Many topical problems in condensed matter physics 
are described by effective Hamiltonians with explicit 
three- or four-body interaction terms, even though the 
underlying Coulomb interaction between particles is only 
a two-body interaction. An example is the Pfafhan state 
of paired composite fermions in the fractional quantum 
Hall effect [ij, which is the exact ground state of a 
three-body interaction this state has been observed 
&t v — I Q. Regular Josephson-junction arrays of a 
p + ip superconductor like Sr2Ru04 can be modeled by a 
multiple-spin lattice Hamiltonian j^] , as can several mod- 
els of frustrated magnetism |^ and superconductivity . 
This paper considers quantum effects on a classical four- 
spin Hamiltonian as an improved model of a p -t- ip su- 
perconducting array, and presents several exact results 
on the resulting 2D quantum phase transition. 

This analysis is based on an exact strong-weak cou- 
pling self-duality for multiple-spin interactions in high di- 
mensions, generalizing the Kramers- Wannier duality 7] 
of the classical Ising model in 2D or the quantum Ising 
chain in ID. This higher-dimensional self-duality contin- 
ues recent developments 0, Q showing that phenomena 
that occur with two-body interactions in one quantum 
dimension, like spin-charge separation, can also be real- 
ized by three- or four-body interactions in two quantum 
dimensions. Superconducting arrays and frustrated mag- 
nets are important examples of this physics because they 
can generate three-, four-, or six-spin interactions with- 
out two-spin interactions, essentially because of unusual 
symmetries. There are many known dualities in d > 2 
that relate a strong-coupling regime of one model to a 
weak-coupling regime of another model, e.g. the dual- 
ity between the 3D Ising model and 3D Z2 lattice gauge 
theory, but self-dualities are quite rare in d > 2 9]. 

We start from the following classical model of Ising 
spins on a square lattice: 
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(1) 



where = ±1, i = 1, ... ,4 are Ising variables at the 
four corners of one face of the lattice, and the sum is 
over all faces. This model was recently introduced 
to understand the effects of frustrating geometric phases 
in a square array of superconducting grains where each 
grain has either p + ip or p — ip order: the state of grain 
i is described by both a phase and an Ising variable 
Si = ±1 that determines the order parameter p + isip. 

The same geometric phases that led to the experimen- 
tal determination of d order in the cuprates [lol lead 
to frustration of the superconductivity unless each pla- 
quette of four grains has an even number of Ising -1-1 
spins (and an even number of —1 spins) 0. The phase 
$ acquired by a Cooper pair moving around a plaquette 
is determined by the states = ±1 of the four grains at 
corners z = 1, . . . , 4: ^ — ^{si + S2 + + S4). This phase 
is equivalent to zero if the plaquette has an even number 
of spins, and otherwise equivalent to tt, which gen- 
erates the local energ y in ijljl proportional to siS2S3S4- 
Josephson weak links 12!| have been made in superfluid 
He"^, and in one phase the symmetries of the weak link 
break the symmetry group down to p ± ip. 

The classical model (QJ also describes the two- 
dimensional "right-angle water" ice model 1131 and maps 
onto to a case of the eight-vertex model P.ll4l|. Note that 
(|T|) is not ordinary 2D Z2 gauge theory, where the sum is 
over bond variables cTi around each face. The overall sym- 
metry group is much smaller for (^: Vol(G") = 2^==+^" 
rather than the full gauge group Vol(G) — 2^^^^ of Z2 
gauge theory. The problem has a one-dimensional 
ground-state degeneracy 2^^+^" even with no physical 
boundary (e.g. on a torus). The model can be solved in 
the thermodynamic limit for all K: its free energy per 
face is just Pf — — log(2coshif), since all the face vari- 
ables can be chosen independently. The model is equiv- 
alent to the ID Ising model and has no phase transition. 

The model becomes quantum-mechanical in the pres- 
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ence of a transverse magnetic field: 



□ 



(2) 



Here again the K interaction is around a plaquette but 
now the spin is a quantum spin-half and the a are Pauli 
matrices. In the superconducting array realization, the 
magnetic field h corresponds to tunneling between the 
two order parameters p ± ip, which will in the limit of 
strong tunneling induce a single real order parameter px- 
Application of pressure in Sr2Ru04 is found experimen- 
tally to drive the system toward a real p state [13 , but 
the explanation of this effect is unclear. The model Q 
is clearly one of the simplest possible 2D lattice quan- 
tum Hamiltonians with four-spin interactions. The main 
weakness of the model (0) is that in the real system at 
low T there is a long-ranged vortex-vortex interaction 
between frustrated plaquettes 0]; the Hamiltonian (0) 
corresponds to treating the core energy of a vortex but 
not its interaction with other vortices, as might be ap- 
propriate at higher temperatures. 

We now show that the quantum model lO has a phase 
transition at zero temperature when K/h is exactly one. 
For simplicity we will give the model's self-duality here in 
its classical form; the same duality can be shown directly 
in the quantum model (0) and interchanges K and h [l6| . 
The 3D anisotropic classical model 
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FIG. 1: The classical anisotropic 3D problem that describes 
the quantum critical point of the model The two types of 
interactions (shaded bonds) are plaquette interactions in the 
planes normal to z, and bond interactions along z. 



]^(cosh Jz + S1S2 suih Jz) 



(7) 



Introduce face variables ku — 0,1 and bond variables 
kh = 0,1, and define cq = coshii', ci = smhK, do = 
cosh Jz , c?i = sinh . Then 



Z = 



{s} ka feb \ □ / 



b h\kh 



(8) 



where □ ranges over all plaquettes in the xy planes, and 
b ranges over all bonds in the z direction, will be shown 
to have a phase transition along the line 



sinh 2/^ sinh 2 Jz = 1. 



(4) 



The connection between coupling constants in the classi- 
cal model and in the quantum model is standard 0| : 



K = aK, 



= tanh(a/i). 



where a is the lattice spacing and M the number of sites 
in the z direction of the classical model, and T is the 
temperature in the quantum case. 

The quantum-classical mapping becomes exact in the 
limits T ^ 0, a ^ 0, and M 00. Knowing the phase 
transition line in the classical model I^J fixes the quan- 
tum transition because in the above limits 

sinh 2k sinh 2J, Ke^'^' = K/h =1. (6) 

The partition function of the above model is 



^ W [(coshX + s^s^s^s^ sinhX) 
{s} □ 




Now the spin sum can be evaluated: for each spin the 
result is 2 if the spin is raised to an even power, and 
otherwise. Z is a constrained sum over the k variables: 



(9) 



Here TV = N^NyN^ is the total number of sites. Each 
site of the original lattice appears via 4 face terms and 2 
bond terms. The constraint is that the sum of the six k 
variables be an even number for every site. 

Now introduce dual variables to solve the constraint. 
The dual spins a are located at the centers of the fun- 
damental cubes of the original cubic lattice. For a site i 
of the original lattice, its four neighboring spacelike faces 
are pierced by four vertical bonds of the dual lattice, and 
for each piercing bond b of the dual lattice fix the rela- 
tion /cq = i(l — crJcTj). Each of the two vertical bonds 
b containing site i pierces a spacelike face □ of the dual 
lattice, and we set ki, = i(l — a^a^a^a^). 

These variables satisfy the constraint since the eight 
dual lattice sites cti , . . . , cts around an original site satisfy 
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+cr50'60'70'8 + trio's + (^2^5 + fsC^T + crA(^%) = 0mod2. 

(10) 

This holds if aU spins are up, and flipping any spin 
changes the sum by an even number. Next we need to 
find how many dual spin configurations correspond to 
one configuration of the k variables. The answer is just 
the size of the gauge group Vol(G") = 2^="+^", since once 
the dual spin configuration is set on a spacelike plane, the 
vertical bonds fix the configuration everywhere else. 
The last step is to calculate the dual couplings. Writing 



Cfc = fc sinh K + {1 — k) cosh 



(11) 



for a face □ of the original lattice, pierced by Ising bond 
b in the dual problem, 

= ^^^4^ cosh k + i^^^ sinh K 

= 4^(1 + CT?(TUanh J*) ^ e^^i^^' (12) 

2 ^ ^ ^ ^2 sinh 2 J* ^ ' 



where tanh J* — e . By the same process 
1 



\/ 2 sinh lk* 



'1 "2 "3 "4 



(13) 



with tanhi^* = e"^"'-. 

Now we combine the above results: the number of 
bonds parallel to z and the number of spacelike faces are 
both N = NxNyNz, so we have up to boundary terms 



z{k,j,) = 



z{k\Jt) 



(14) 



Vol(G") sinh(2 J*)^/2 sinh(2ii:*)^/2 
More precisely, we have shown that in the limit iV — )■ oo, 



iogZ(x,j,) \ogZ{k\j:) 1 



iV 



N 

logsinh(2j; 



logsinh(2ii:*) 



The duality relation for the couplings is 



sinh(2J*) 



sinh(2if*) = ^ 



sinh(2X) 



sinh(2J2)' 



(15) 



(16) 



This determines the entire phase boundary line in the 
(i^, Jz) plane. It also shows that there is a hidden sym- 
metry between the couplings K and Jz- The same con- 
struction leads to an exact duality of the 4D model with 
eight-spin interactions around spacelike cubes, plus a 
bond interaction in the fourth dimension; and similarly 
for any d > 2. The self-dual point is i^T = = Kc = 
i log(l -I- V2). If there is only one phase transition in the 
model at finite coupling, it must be at the self-dual line. 
The self-duality of this model is similar to the classical 
anisotropic self-dual Villain Zn models studied in 0. 

We remark that the model can be solved if the sys- 
tem has only one row of spacelike plaquettes: the bond 



variables = Sx.iSx,2 become spins in an anisotropic 2D 
Ising model, and with e^"^ = cosh 2 J^, 



logZ 



{K,Jz 



logZ],%^^{2k,J') 
NxNz 



(17) 



If there is a single second-order transition in the 3D 
classical model, then there is a second-order transition 
'At K — h 'm the quantum model. The spontaneous or- 
der that develops across the transition is unconventional 
"bond order" . At /i = 0, the ground states all have long- 
range order along every row of the vertical bond erf (Tf_|_y , 
and of horizontal bonds along each column. Another way 
to describe the h — Q ordered state is using the Wilson 
operator for the product of spins around a loop £, 



(18) 



In any ground state, {W{C)) = 1. We will use bond 
order below to make a connection between the order pa- 
rameter and the excitation spectrum. The two orders are 
connected since a closed loop contains an even number 
of bonds from each row and column. 

A Peierls-type argument can be used to show that there 
is an ordered phase of the classical 3D model at low tem- 
perature, and hence at least one phase transition. We 
have performed Monte Carlo and high-temperature se- 
ries calculations to check whether the above model has a 
single second-order phase transition (which must then lie 
at the self-dual point Kc). The results are consistent with 
this picture, but the single-spin Monte Carlo algorithm 
becomes very slow close to criticality, as in 0|. 

The high-temperature series proves that this model 
does not have the same free energy per site as an Ising 
model, as might have been suspected since the classical 
2D model with face interactions ^ has the same free 
energy per site as the ID two-spin Ising model. The first 
terms for the classical symmetric model {K — J^) are 



" dT N{dTf 



(19) 

This differs from the d>2 two-spin Ising model at order 
k'^ and the d=l two-spin Ising model at order k^ . 

The quantum model ^ has an infinite but non- 
extensive set of conservation laws: along any one of the 
Ny rows, say 7?., the product 



O 



(20) 



commutes with the Hamiltonian and similarly for 
each of columns. These N^+Ny conserved quantities 
are related to the ground-state degeneracy in the ordered 
state. There are 2^=^+^" sectors of the theory, labeled 
by the eigenvalues of the operators H20|) . In the large- 
h phase, there is a single ground state invariant under 
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Nx + Ny transformations that each act on all the spins 
in one row or column via 

Cx CTx, (Jy — CTy, CTz —<Jz- (21) 

At the transition K ~ h va. the thermodynamic limit, 
these symmetries are spontaneously broken and there are 
2N^+Ny degenerate ground states once K > h. In a fi- 
nite system, there is no spontaneous symmetry break- 
ing and ground-state degeneracy on the ordered side 
K > h will be split by an amount exponentially small 
in -nim.{Nx,Ny). The breaking of these many symmetries 
at a single transition occurs because of the infinite num- 
ber of conservation laws. In the language of hard-core 
bosons, the charge along each row or column is conserved 
modulo 2, just as in the Bose metal model discussed in 0. 

Now we consider the excitation spectrum in the two 
phases of the quantum model. Recall the familiar quan- 
tum Ising chain (d = 1) the lowest excitation in 
the ordered phase with periodic boundary conditions is 
to flip one spin from the ground state, so the first excited 
state contains two bad bonds. These two bad bonds can 
be separated into two kinks, with a string of flipped spins 
between them. In the large h limit, all the spins point 
along X and the one-particle state is a flipped spin with 
momentum and kinetic energy, with 

ek= Kh [2-2/ hcos{k)+0{l/h^)\ . (22) 

For our model in 2D, in the large K limit the system 
stays near the ground state manifold, and the first excited 
state locally connected to the ground state is obtained by 
flipping one spin. This flipped spin results in four bad 
plaquettes (or four Z2 vortices). As in the d — 1 case, 
this excitation can disintegrate into four fractional exci- 
tations, which become four vertices of a rectangle with 
all the interior spins flipped. 

The large-/i limit requires more attention. In the x ba- 
sis, it is clear that the single flipped spin is nondispersive 
because of the Z2 conservation laws. The single flipped 
spin in the large-/i limit is equivalent via duality to a 
single bad plaquette in the large-if limit. Instead, the 
lowest mobile state is a flipped bond: the flipped bond 
can hop in one direction, the direction perpendicular to 
the bond, and its dispersion relation is the same as that 
in the d — \ case (|22ll . These one-particle states can scat- 
ter off each other, and because they hop unidirectionally, 
this kind of scattering is very similar to ID scattering. 
The ordered state at small his a condensation of bonds. 

Direct experimental observation of the order- 
parameter state of a superconducting grain is currently 
only possible for isolated superconducting grains. How- 
ever, there are several experimental signatures of the 
bond ordering predicted above, even though in the 
ground state the system has no frustrating fluxes and 
hence is uniform from the superconducting point of view. 
The gapless fluctuations near the transition will modify 



the specific heat and transport in the system: transport 
will be attenuated near the critical point by scattering 
off the fluctuations, while specific heat will show a 
peak. The true spin action 'm a, p + ip superconducting 
array is more complicated and includes long-ranged 
interactions between frustrated plaquettes 0] , but it has 
the same gauge symmetries as (0). If grain orders could 
be measured at the array boundary, the directionality 
of the predicted bond ordering would be seen, and the 
Ising variables of each neighboring pair of grains develop 
a symmetry-breaking correlation at the transition. 

An unresolved question is whether ((^jl is a free-fermion 
model like ID quantum Ising. We have not found addi- 
tional conservation laws beyond ordinary symmetries and 
the infinitely many Z2 laws described above. A solution 
of this 2D quantum model would give further insight into 
the physics of explicit four-point interactions. 
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